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ABSTRACT 

We present a physically motivated semi-analytic model to understand the clustering 
of high redshift Lyman break galaxies (LBGs). We show the model parameters con- 
strained by the observed luminosity function of high-z LBGs can be used to predict 
large scale (i.e 9 ^ 80 arcsec) bias and angular correlation function of galaxies. These 
predictions are shown to reproduce the observations remarkably well. We then adopt 
C"| . these model parameters to calculate the halo occupation distribution (HOD) using 

O [ the conditional mass function. The halo model using this HOD is shown to provide 

a reasonably good fit to the observed clustering of LBGs at both large (9 > 80") 
and small (9 < 10") angular scales for the whole range of z = 3 — 5 and limiting 
magnitudes. However, our models underpredict the clustering at intermediate angular 
scales, where quasi-linear effects are important. The average mass of halos contribut- 
ing to the observed clustering is found to be 5.7 x 10 11 M Q and the characteristic 
mass of a parent halo hosting satellite galaxies is 10 12 M Q for a limiting absolute 
magnitude of —20.5 at z = 4. For a given threshold luminosity these masses increase 
with increasing z and at any given z these are found to increase with increasing value 
of threshold luminosity. We find that the average number of satellites in a parent halo 
is smaller than 1. Thus the small scale clustering detected at high-z is most likely 
to be dominated by LBG pairs rather than few large clusters. Large scale clustering 
predictions have only weak dependence on our model parameters and we may need 
measurements with better than 10% accuracy to be able to constrain the cosmological 
parameters. However, we find that the small angular scale clustering is quite sensitive 
to changes in the duration of star formation in a halo and hence could provide a probe 
of this quantity. The present data favor star formation in a halo lasting typically for 
a dynamical time-scale. 

Key words: cosmology: theory - cosmology: large-scale structure of universe - galax- 
ies: formation - galaxies: high-redshift - galaxies: luminosity function - galaxies: statis- 
tics - galaxy: haloes 



1 INTRODUCTION 

Over the past decade there has been a growing wealth of 
observations probing the properties of high redshift galax- 
ies. Various surve ys, using the Lyman break col or se 



les. various surve ys, using tne .Lyman preak color se- 
lection technique (IMadau et af] 1 19961: ISteidel et al] 1 19961 : 



lAdelberger et al.lll998l : ISteidel et al.lll998l ). have detected a 
substantial number of high redshift galaxies, up to z ~ 10. 
This has resulted in reasonably good estimates of luminos- 
ity functions (LF) of these Lyman break galaxies (LBG) 
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up to z 



( Bouwens et ai1l2007l. 120081 : iReddv et al1l2008l : 



201C 



2 ~ 5 (lOuchi et al. 


20041. 20051: Kashikawa et all 


2006 


Lee et alj 2006: Hildebrandt et al.ll2009l: ISavov et all 


2011 


Bielbv et al.l 201lf). It is important to explain these obser- 



vations and understand their implications for galaxy forma- 
tion. 

In the hierarchical model of structure formation galax- 
ies form in virialized dark matter halos. These inturn result 
from the growth and gravitational collapse of initial Gaus- 
sian density perturbations. Thus the statistical properties of 
galaxies are determined by the statistics of the parent halo 
population, given a model for how they form inside these ha- 
los. The properties of dark matter halos are quite well under- 
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stood using N-body simulations jSpringel et al.ll2005l ) and 
analy tical models like the h alo model of large scale struc- 
ture (|Coorav fc Shethll2002l ). These approaches provide the 
abundance, spatial distribution and merger history of dark 
matter halos. Numerical simulations also s uggest a possi- 
ble u niversal dark matter halo density profile (|Navarro et al.l 
Il997h . Given the above inputs on dark matter halo proper- 
ties and a specific model of galaxy formation inside these 
halos, it is possible to explain the two major observables of 
galaxies, their luminosity function and clustering. In addi- 
tion, such models can throw light on the complex physics of 
galaxy formation, such as rate and duration of star forma- 
tion, feedback mechanisms etc. 

There has been extensive modelling of the luminos- 
ity function s and clustering of gala x ies at low redshift s 
|Yang et al.l 120031 ; IZheng et all |2009| ; IZehavi et all l201lf ). 
Several studies on understanding the LF of high red- 
shift LBGs have also been carried out (|Benson et al.l [20031 ; 
IStark et al.ll2007l ; iKhochfar et all 120071 '). We have also ear- 
lier extensively studied simple physical m odels of galaxy for- 
matio n to understand the LFs of LBGs (|Samui et al. I l2007l . 
l2009h . These studies have reasonably reproduced the LFs 
of high z LBGs from z — 3 — 7. In addition they constrain 
the efficiency of star formation and its duration in LBGs, 
and have al so been used to set tight limits on the neu- 
trino mass (jjose et al.l l201lh . We now wish to examine if 
our simple physical model for galaxy formation, combined 
with the halo model, can also explain the clustering of the 
high z LBGs. Recent work on understanding LBG cluster- 
ing is base d on empirical or parametrized models of galaxy 
formation jBullock et al.ll2002al; lHamana et all 1200412009 ; 



IConrov et al.ll200d : lLee et al.ll2009l )7 In our work we explore 
this issue using a physically motivated model of star forma- 
tion inside dark matter halos. 

In order to model the clustering of LBGs, we 
choose an approach similar to the convent ional halo 
occupation distribution (HOD) metho d lljing et al 



occupation distributio n (fiUD) metno d IIJmg et al 
19981; ISeliakl |2000|; IScoccimarro etail l200ll; iBuIlock et a! 



2002bl; iBerlind fc Weinberg! |2002|; IBerlind et all 120031: 



Kraytsov et al.ll2004l : IZehavi et al.ll2004l . 120051 : IZheng et ail 
20051 ). In our model, the luminosity of any galaxy is 



obtained from a physical model of star formation rate 
which depends on the mass and age of the hosting halo. 
To begin with we assume that each halo can host at most 
one visible galaxy. Under this assumption we show that our 
prescription of star formation can simultaneously explain 
both the luminosity function and large scale clustering 
of high redshift LBGs. Further we extend our approach 
using a physically motivated HOD prescription; where 
satellite galaxies harbored inside dark matter halos are 
also incorporated. It should be noted that, unlike earlier 
works on high-z LBGs, we separate the HOD into the 
sum of central and satellite galaxies, in acco rdance with 
the previou s models in low red shift universe (|Zheng et alj 
120051 . 120091 ; IZehavi et all l201lh . To calculate the number 
and luminosity of the sat ellite galaxies we emp loy the 
conditional mass function (|Coorav fc Shethl 12002} ) . Using 
our approach, we find that one can explain both the UV 
LFs and luminosity dependent clustering of LBGs and 
obtain useful physical insights into galaxy formation. 

The organization of this paper is a follows. In the next 
section we describe our physical model for computation of 



the LF of LBGs. In Section 3 we focus on the clustering of 
LBGs on large angular scales. We then turn to our physi- 
cally motivated model to calculate the central and satellite 
contributions to the HOD and use these to obtain total an- 
gular correlation functions at all angular scales. Section 5 
presents a comprehensive comparison of the total angular 
correlation function computed in various models with ob- 
servations. A discussion of our results and conclusions are 
presented in the final section. For all calculations we adopt 
a flat ACDM universe with cosmological parameters consis- 
tent with 7 year Wilkins on Microwav e Ani sotropy Probe 
(WMAP7) observations (|Larson et al.l l201ll) . Accordingly 
we assume fi m = 0.27, = 0.73, fi b = 0.045, h = 0.71, 
n s — 0.963 and as = 0.801/i _1 Mpc. Here Qi is the back- 
ground density of any species 'i' in units of critical density 
pcric- The Hubble constant is Ho = WOh km s _1 Mpc -1 



2 THE STAR FORMATION RATE AND 
LUMINOSITY FUNCTION 

Our aim is to construct a self consistent semi-analytical 
model, that can explain the luminosity functions and clus- 
tering of high redshift LBGs. In this section, we briefly 
recall the semi analytical treatment of ISamuT et al.l (|2007l ) 
(hereafter SSS07) to model luminosity fu nctions of high red- 
shift LBGs [See a lso ISamui et all (|2009l ) (hereafter SSS09); 
IJose et all (|201ll )] before presenting our prescription to 
calculate angular correlation functions at large and small 
scales. 

In the models of SSS07, the star formation rate (Msf) 
in a dark matter halo of mass M collaps ed at redshift z c 
and observed at redshift z is given by (see. lChiu fc Ostrikerl 
120001 ; IChoudhurv fc Srianandll^Ool ) , 



M S f(M,z,z c ) 



t(z)-t{z c ) 



(i) 



x exp 



t(z)-t(z c ) 



where, /* is the fraction of the total baryonic mass that is 
converted into stars over the entire lifetime of the galaxy. 
Here, td y n(z c ), is the dynamical time scale of a halo collaps- 
ing at z c (Eq. (3) of SSS07) and k is a parameter which 
governs the duration of the star formation activity. Further, 
t(z) is the age of the universe at redshift z; thus t(z) — t(z c ) 
is the age of the galaxy at z. 

The stars are formed with a normal Salpeter IMF in 
the mass range 1 — 100 Mq. The pop ulation synthesis code 
STARBURST99 (|Leitherer et al.lll999l ) is used to obtain the 
rest frame luminosity (Z1500) at 1500 A as a function of time 
of a galaxy undergoing a burst of star formation. The as- 
sumed star formation rate of a galaxy, as given in Eq. (flj, 
is then convolved with this burst luminosity to get the time 
evolution of the luminosity, L1500, of an individual star form- 
ing galaxy (See Eq. (6) and Figure 1 of SSS07) 



(2) 



£i50o(i) — I Msf(T — t)1 15 oo(t)cIt. 



Due to dust absorption, only a fraction (1/t/) of L1500 
produced by the stars manages to escape the galaxy. This 
luminosity (L = L1500/77) is then converted to a standard 
absolute AB magnitude Mab, using the equation given 
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Figure 1. Comparison of obse rved UV LF o f LBG s at three diffe rent redshifts wi t h our best fitted model predictions. The observed data 
points and error bars are from iReddv et all IMPOST ) (for z=3) and lBouwens et al.l J2007T) (for z=4 k, 5). 



by (|Oke fc GunrJ Il983p . to enable direct comparison with 
the observed data. Having obtained the Mab of individual 
galaxies we can compute the luminosity function $(Mab, z) 
at any redshift z using, 



$>(M AB ,z)dM AB 



dn(M(M A 



dM 



dhx 



dM A i 



dz c 



dM AB . 



dL 15 



(3) 



Here dn(M, z c )/dz c = n(M, z c )dt/dz c , and h(M,z c )dM is 
the formation rate of halos in the mass range (M, M+dM) at 
redshift z c . S SS09 modelled th i s form ation rate as the time 
derivative of ISheth fc Tormenl |l999h (hereafter ST) mass 
function as they are found to be the good in reproducing the 
observed LF of high-z LBGs. Therefore we use n(M, z c ) = 
dnsr{M, z c )/dt where nsr(M, z c ) is the ST mass function 
at z c . Also note that we use the notation n(M) for dn/dAI 
for convenience. 

Star formation in a given halo also depends on the cool- 
ing efficiency of the gas and various other feedback processes. 
We assume that gas in halos with virial temperatures (T v i r ) 
in excess of 10 4 K can cool (due to recombination line cooling 
from hydrogen and helium) and collapse to form stars. How- 
ever the ionization of the IGM by UV photons increases the 
temperature of the gas thereby increasing the Jean's mass 
for collapse. Thus in ionized regions, we incorporate this 
feedback by a complete suppression of galaxy formation in 
halos with circular velocity v c ^ 35 km s - 1 and no sup 
sion with v c ^ V u = 95 km s _1 (jBromm 
intermediate circular velocities, a linear fit from 1 to i s 
adopte d as the suppression f actor iBromm fc Loebl (|2002l) ; 
see also lBenson et al.l l|2002l ): lDiikstra et al.l (|2004h . SSS07J. 
SSS07 found that this feedback mechanism naturally leads 
to the observed flattening of the LF at the low luminosity 
end. 

In our models, we also incorporate the possible Active 
Galactic Nuclei (AGN) feedback that suppresses star for- 
mation in the high mass, by multiplying the star formation 



apres- 
Loebll2002h . For 



rate by a factor [1 + (M/Magn)^]^ 1 . This decreases the star 
formation activity in high mass halos above a characteristic 
mass scale M a0 n, which is believe d to be ~ 10 12 Mq (see 
iBower et aHl2006l ; iBest et al.ll2006h . In our models we con- 
sider P = 0.5 and M agn as a free parameter. Note, SSS07 
used AGN feedback with j3 — 3, but suppressed the halo 
formation rate instead of star formation rate as we do here. 

A crucial parameter of our model is f*/rj which gov- 
erns the mass to light ratio of the galaxies at any given 
redshift. This is fixed by fitting the observed luminosity 
function of LBGs using \ 2 minimization. In this way our 
physically motivated model for star formation gives the re- 
lationship between the halo mass (M) and the luminosity 
(L) of the galaxy it hosts. It is important to note that our 
prescription of star formation naturally introduces a scatter 
in the M-L relationship because halos of mass M forming 
at different redshifts, produce different luminosities at the 
redshift of observation. SSS09 showed that this M-L rela- 
tionship can successfully explain the luminosity functions of 
high redshift LBGs. In FigQ] we show the observed lumi- 
nosity function of LBGs together with our best fitted model 
results for three different redshifts. As a fiducial model we 
have chosen k = 1.0 at all redshifts and three different values 
for Magn, 1.0 x 1O 12 M , 1.8 x lO 12 Af , and 2.0 x 1O 12 M at 
redshifts 3, 4 and 5 respectively. With these parameters the 
observed luminosity functions is well reproduced for /» /rj of 
0.038, 0.037 and 0.032 respectively for z = 3, 4 and 5. The 
X 2 corresponding to the best fit luminosity function at red- 
shifts 3, 4 and 5 are 6.09, 8.3 and 11.67 (with corresponding 
reduced \ 2 of 0.87, 0.83 and 1.66) respectively. 



3 CLUSTERING OF LBGS AT LARGE 
ANGULAR SCALES 

We couple the semi-analytic models described above with 
the halo model to compute the correlation function of high- 
z LBGs. In order to calculate the galaxy-galaxy correlation 
function on all scales one requires a full knowledge of the 
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halo occupation distribution (HOD), which describes the 
conditional probability P(N\M) for N galaxies of a given 
type to reside inside a halo of mass M (|Bullock et al.ll2002bl ; 
iBerlind fc Weinberd[20o3 ). On scales much bigger than the 
virial radius of a typical halo, the clustering amplitude is 
dominated by correlation between galaxies inside separate 
halos. On the other hand, on scales smaller than the typical 
virial radius of a dark matter halo, the major contribution 
to galaxy clustering is from galaxies residing in the same 
halo. These separate contributions to two point correlation 
function are called 2-halo and 1-halo terms respectively. In 
this section we concentrate on the 2-halo term. 

The first moment of halo occupation or the mean num- 
ber of galaxies of a given type inside a parent halo has two 
separate contribut i ons from central a nd satellite galaxies 
|Zheng et al.ll2005l ; IZehavi et al.l 120111 ). In the calculation 
of LF in the previous section we have assumed that each 
halo hosts a single star forming galaxy. However, the de- 
tectability of this galaxy depends on its age, the limiting 
luminosity of the observations and feedback processes intro- 
duced in the previous section. Neglecting one halo term is a 
good approximation for correlation functions on large scales 
because (i) on largest scales the clustering is insensitive to 
the galaxy distribution inside a halo and is dominated by the 
two halo term (ii) using our model, we later show that, for 
the LBGs we consider, the total number of satellite galax- 
ies of a particular luminosity is much less than the total 
number of central galaxies of the same luminosity in a large 
volume. Thus on large scales the contribution to clustering 
due to satellite galaxies is not significant compared to the 
clustering of central galaxies. 

The two halo galaxy power spectru m on large scales 
(sma ll k), assuming linear bias, is given by (|Coorav fc Sheta 
l2002h 



P g 2h (k,z) = b 2 g (k,z)P lin (k,z), 



(4) 



where Pu„(k) is the linear dark matter power spectrum. 
The scale dependent, galaxy number weighted, halo bias or 
in short galaxy bias, b g (k,z), is given by 



b g (k,z) = 



n g (z) 



b(M,z)n(M,z)u(k,M,z)dM. (5) 



As before, n(M,z) is the ST halo mass function, n g (z) = 
f n(M, z)dM is the number density of galaxies and b(M, z) 
is the mass dependent halo bias fac tor provided by the fit- 
ting function of Sheth and Tormen (jSheth fc Tormenlll999l : 
ICoorav fc Shethll2002l ). It has the following functional form 



b(M,z) 



qu(M,z) - 1 
5Jz) 



+ 



2p/5 c (z) 



l+(qu(M,z))p-' 



(6) 



where, v(M,z) — 8 c (z) /o(M), a(M) is the linearly extrap- 
olated rms density fluctuation on any mass scale M and 
& c (z) is the critical density required for collapse at z. Here 
5c\z) = D(z)S c {z = 0) with S a (z = 0) = 1.686. Also 
we use p = 0.3 and q = 0.7 07 as given by Sheth and 
Tormen (|Coorav fc Shetbll2002r i. In Eq. © u(k,M) is the 
Fourier transform of dark matter density profile normalized 
by it's mass, ie u(k,M,z) = p(k, M, z) / M . We assume the 
Navarro, Frenk and Wh ite (NFW) profile (jNavarro et alj 
1 19971 : iBullock et al.llioOlh for the dark matter density dis- 
tribution in halos. On scales much greater than the virial 
radius of a typical halo of mass M, u(k,M,z) = 1. Thus 



galaxy bias, b g (k,z), on large scales (much larger than the 
virial radius of typical collapsed halos) or small k, will be 
independent of k. 

The actual clustering observations are for a galaxy sam- 
ple with an apparent magnitude below some threshold (or 
a corresponding lower luminosity threshold Lth)- So b g (k, z) 
in Eq. ((Sj should be calculated for only galaxies which have 
luminosity greater than L t n- We use our model described in 
Section [2] for calculating the time dependent luminosity of 
a galaxy hosted by a dark matter halo of mass M. This can 
then be used to obtain the galaxy bias for LBGs with the 
luminosity threshold L t h, at any redshift of observation z 



b g (k,Lth,z) 



n g (L th ,z) J 



dMb(M, z)u{k, M, z) x 



dz c Q{L{M,z c ,z) - L th ) 



dn(M, z c ) 
dz c 

(7) 

Here, L(M, z c , z) is the luminosity of a galaxy of mass M at 
z, formed at z c . The theta function, Q(L(M, z c ,z) — L t h), 
in the above equation ensures that the galaxy bias will have 
contribution only from those galaxies formed at z c and that 
shine above the threshold luminosity L t h at z. If we do not 
have this constraint imposed by the O function, the inte- 
gral over z c is just f°° dz c dn(M, z c )/dz c = n(M, z) and the 
Eq. (0 reduces back to Eq. ([SJ. Also n g (z, Lth) is the num- 
ber density of galaxies having luminosity in excess of the 
limiting luminosity, L t h, and is given by, 



n g (L t h,z) 



dM 



dz c Q(L(M, z c ,z) - Lth) 

dn(M, z c ) 
dz c 



(8) 



In our prescription of star formation, at any given red- 
shift z, a galaxy hosted by a halo of mass M will shine with 
the same luminosity when its star formation is either in the 
rising phase or in the declining phase. As a result there are 
two redshifts of formation z\(L t h, M, z) and ziiLth-, M, z), 
such that the galaxy will produce an observed luminosity 
Lth at z. We demonstrate this in Fig. [2l where we have plot- 
ted the halo mass M that can host a galaxy of luminosity 
L t h (corresponding to an apparent magnitude m = 25) at 
2 = 4 against the redshift of formation of the halo. The figure 
clearly shows that halos of any mass (above a minimum mass 
M m in) formed at two different redshifts can shine with the 
same luminosity at z = 4. For example, when M = 1O 12 M0 
these two redshifts are approximately 4.24 and 6.01 (shown 
by the abscissa of intersection of the horizontal red dotted 
line with the curve). A halo of mass 1O 12 M formed be- 
tween these redshifts will shine with a magnitude brighter 
than Mab = —20.1 at z = 4. For the case illustrated above, 
this minimum mass is roughly 4.8 x 10 11 Mq and this halo 
has to be formed at z c ~ 4.8 for it to have a luminosity L t h 
at z = 4. In our models, the exact values of zi , Z2 and M m in 
for a given observed luminosity will depend on f*/rj that we 
fix by fitting the luminosity function as discussed before. 

Our model of star formation ensures that halos of 
mass M > M min collapsing between z\ (Lth , M, z) and 
z%[Lth, M, z) will always shine brighter than L t h at z, while 
galaxies formed outside these intervals do not. Thus the O 
function in Eq. [7j is unity for M > M m i n and zi ^ z c ^ Z2 
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Figure 2. Masses of the dark matter halos that shine with an ap- 
parent magnitude m = 25 (or absolute magnitude Mab = —21.1) 
at z = 4 as a function of their formation (collapse) redshift. We 
can clearly see that halos of a particular mass formed at two 
different redshifts can shine with the same brightness at z = 4. 
We refer these two redshifts as zi and 22. For example, when 
M = 1O 12 M0 (shown by dotted red horizontal line ) these two 
redshifts are approximately 4.24 and 6.01 The figure also shows 
the minimum mass of the galaxy that can produce an apparent 
magnitude 25 at z = 4 (in red horizontal solid line). This mass is 
roughly 4.8 X 10 11 Mq and is formed at z c ~ 4.8. 



and zero otherwise. Therefore b g (k, L t h 
written as, 



2) in Eq. ([7]) can be 



b g (k, 



1 



n g {L th , z) 



dMb(M, z)u(k, M, z) 



£ 



dz, 



dn(M, 



dZr 



(9) 



At this point, given M min ,z 1 and Z2 one can compute 
b g (k, Lth, z). As these parameters can be fixed by the param- 
eters governing the star formation in a halo, the observed 
luminosity function alone in principle will allow us to predict 
b g (k, Lth, z) uniquely. 

We present our results for the fiducial set of model pa- 
rameters that reproduces the observed luminosity function. 
In Table [5] we present the predicted asymptotic (k — ► 0) 
value of b g (Lth, z) at three different redshifts (and three lu- 
minosity ranges each) using f*/r] values obtained by fitting 
the observed luminosity function of LBGs at these redshifts. 
This table also gives apparent magnitude cut-off (m), cor- 
responding absolute magnitude cut-off (Mab), M m in and 
n g (L t h,z) as given in Eq. 8. 

We compute the luminosity dependent galaxy power 
spectrum, P g h (k, L t h, z), by substituting b g (k, L t h, z) from 
Eq. © into Eq. @. We then have 

Pg h (k,L t h,z) = b 2 g (k, Lth, z)Pu n (k, z). (10) 

The corresponding luminosity dependent two point correla- 
tion function of galaxies with lum inosity greate r than Lth 
at z, can now be calculated using (|PeeblesHl980T ) 



i 2 g h {r,z,L t h) = 



dk 2 sin(fcr) 2h . . . . 



Table 1. Asymptotic values of b g (L t h, z) predicted uniquely by 
using model parameters that best fit the high z LFs. These pa- 
rameters are (i)/*/»j, an indicator of the light to mass ratio at any 
redshift (ii) M m ; n {Lth , z ), the minimum mass of a galaxy that can 
shine brighter than a given luminosity threshold L t h at redshift 
of observation, (hi) n g (z,L t h), the number density of galaxies 
with luminosity greater than L t h in units of 10 — 4 (/i/Mpc) 3 and 
bg(L t h,z), the luminosity dependent galaxy bias are tabulated 
for three limiting apparent magnit udes at each redshift. The last 
column is the galaxy bias given bv lHildebrandt et al.l d2009h . af- 
ter correcting for the larger erg adopted by them. The parameters 
obtained here are used to compute angular correlation function 
on large scales. 



z 


f*/v 


m 


M AB 


M mm /M Q 


n g 


b 9 


b g (H) 






24.5 


-21.1 


6.0 x 10 11 


9.37 


3.46 


4.50 


3 


0.038 


25.0 


-20.6 


3.8 x 10 11 


18.0 


3.16 


3.22 






25.5 


-20.1 


2.4 x 10 11 


33.4 


2.90 


2.67 






25.0 


-21.1 


4.8 x 10 11 


4.89 


4.56 


5.14 


4 


0.037 


25.5 


-20.6 


2.9 x 10 11 


10.9 


4.16 


4.25 






26.0 


-20.1 


1.8 x 10 11 


23.1 


3.81 


3.51 






25.5 


-21.0 


4.3 x 10 11 


1.40 


6.03 


7.81 


5 


0.032 


26.0 


-20.5 


2.7 x 10 11 


3.67 


5.46 


6.08 






26.5 


-20.0 


1.7 x 10 11 


9.01 


4.96 


5.09 



2tt 2 



We compute luminosity dependent angular correlation func- 
tion w(9, z) fr om the spatia l correlation function using Lim- 
ber equation \ Peebles! 1 1980h 

w(6,z)= I dz N(z) I dz" N(z")Zg(z,r(6;z',z")) 
Jo Jo 

(12) 

where r(9;z',z") is the comoving separation between two 
points at z and z" subtending an angle 6 with respect to 
an observer today. Here we have also incorporated the nor- 
malized redshift selection function, N(z), of the observed 
population of galaxies. In Eq. (|12|l we neglect the redshift 
evolution of clustering of the galaxies detected around z. 
Hence the spatial two-point correlation function £ 9 (r, z) is 
always evaluated at the observed redshift. 

In Fig. [3] we over plot the angular correlation func- 
tion computed in this way using the predicted value of b g 
from our models that produce best fits to the LF, on the 
observed angular t wo po int correlation functions given by 
iHildebrandt et all (|2009l ). We have also plotted the angu- 
lar correlation function of dark matter density (obtained 
by putting b g (k,L t h,z) = 1) in dotted-blue curves. To 
compute our corr e lation function we used the N(z) from 
IHildebrandt et al.l (|2009l ). which are kindly provided by the 
authors (BCsim redshift distribution; see Table 4 and Fig- 
ure 5 of the paper). It is clear that at large angular scales 
(i.e 6 ^ 80 arc sec), where the linear approximation used 
here is valid, our model predictions match well with the ob- 
servations at different redshifts and different luminosity (or 
apparent magnitude) thresholds. 

The large scale galaxy bias is determined by 
IHildebrandt et all (|2009l ) from their data by making power 
law fits to observed galaxy correlation functions, and com- 
paring the corresponding galaxy variance at 8/i _1 Mpc 
(c"8,g) w hh as computed from the dark matter power spec- 
trum. In the last column (column 7) of Table. [2] we show 
the large scale bias b g (H) of LBGs obtained this way by 
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Figure 3. The large scale angular correlation function of LBGs at redshifts 3, 4 and 5 for various limiting magnitudes. Each row 
corresponds to a particular redshift, which is labelled in the first panel of that row. In each row there are three panels showing the 
clustering predictions for galaxies with three limiting magnitudes that are labelled in each panel. The blue solid curves are our model 
predictions of galaxy angular correlation functions. Also shown in blue do tted lines are the angular correlation functions of dark matter 
densi ty. The data points and error bars shown by black triangles are from Hildcbran dt et al. and those shown by red squares are 

from lOuchi et al] feoOEft . 



IHildebrandt et all |2009h . after correcting for the larger ag 
adopted by them. The large scale bias b g that we pre- 
dict from our models ag rees well with b g (H) estimated by 
IHildebrandt et all (|200a i from their power law fit, for all but 
the brightest of the LBG samples. For the brightest sample 
of LBGs at all redshifts the bias we predi ct is systematically 
lower. Nevertheless, in agreement with IHildebrandt et all 
l|2009l ). we also find that (i) at a given z the bias increases 
with increasing luminosity and (ii) for a given L t h the bias 
increases with increasing z. However at any given z the 
spread in b g as a function of Lth i s less in our case com- 
pared to that of IHildebrandt et all (120091 ). 

We also see from Fig. [3] that the two halo term is not 
able to account for the strength of clustering seen on small 
angular scales given by 9 < 50" at redshift 3 and 9 < 80" 
at redshifts 4 and 5. However we have not yet included the 
1-halo term in the correlation function, arising from correla- 



tion between galaxies within any given halo. We will address 
this issue in the Section [4] 

3.1 Sensitivity of large angular scale correlation 
function to change in parameters 

Here we present the sensitivity of the predicted large scale 
correlation functions to changes in astrophysical and cosmo- 
logical parameters. To do this, we considered models, where 
a single astrophysical or cosmological parameter is assigned 
values around it's fiducial value, keeping all other parame- 
ters fixed to their fiducial value. 

3.1.1 Sensitivity to Astrophysics 

We begin by showing in Fig. [4] sensitivity of large cluster- 
ing predictions to the assumed values of, k and M agn - Note 
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Figure 4. The percentage change is angular correlation function 
when one vary n (top panel) and M a gn (bottom panel) around 
the fiducial model at z = 3. Here the plots are for threshold 
apparent magnitude 25. For each value of k and M a an we use 
the appropriate /* /n values that best fit the observed luminosity 
function. 
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that, in each case the f*/rj gets automatically fixed when we 
fit the observed luminosity function. The results are shown 
at z = 3 and m = 25. Any change in astrophysical parame- 
ters will not alter the amplitude and shape of dark matter 
correlation function. However in principle they can change 
the luminosity dependent galaxy bias by simply changing 
the luminosity of a galaxy. Hence on large scales the astro- 
physical parameters will affect clustering predictions only 
through galaxy bias. 

As we can see explicitly in this figure, the changes in 
astrophysical parameters produce negligible effects on the 
clustering on large scales. This is expected for M agn be- 
cause of the following reason. For z — 3 we have chosen 
the parameter M agn ~ 1 x 10 12 Mq so that the predicted 
LF matches with the observed data. In the plot we varied 
Mag„ between 0.7 ^ M ag n/W 12 M Q s: 1.3. From Table[2]it 
is clear that, M m in, the minimum mass cutoff in the inte- 
gral to calculate large scale bias (see Eq. [9]) is always much 
less than 10 12 Mq. Since the number density of halos decays 
exponentially with a 2 (M) at these mass scales, the major 
contribution to b g (z,Lth) in Eq. [9] comes from mass scales 
much lower than M agn - Hence we can conclude that the 
value of Magn used in our models will not be that sensitive 
to clustering on large scales. 

From Fig 0] it is also clear that the change in re is not 
affecting the predicted clustering at large scales. We know 
that a change in re is not affecting the number density of 
halos of a particular mass. However for a fixed f*/rj a small 
re shortens the duration of star formation activity resulting 
in an increase in the luminosity of galaxies. This will shift 
the predicted total LF more or less along the luminosity 
axis. While fitting the observed luminosity function this shift 
with respect to our fiducial model is nullified by changing 
ft/rj. This means that even when the re is decreased and the 
star formation rate is enhanced the luminosity of a galaxy 
hosted by a halo of a given mass is almost unchanged, and 



Figure 5. The percentage change is angular correlation function 
when one varies the cosmological parameters around the fiducial 
model at z = 3. Here the plots are for threshold apparent mag- 
nitude 25. In each case, we use the appropriate ft/17 values that 
consistently reproduce the observed luminosity function. 



hence the luminosity dependent large scale clustering is not 
significantly affected. 

3.1.2 Sensitivity to Cosmology 

On large scales values of cosmological parameters can affect 
two point correlation function by changing the amplitude 
and shape of dark matter power spectrum and also by mod- 
ifying the halo mass function and halo bias and thereby al- 
tering the galaxy bias. To see these effects more clearly, we 
have plotted in Fig.[S]the percentage differences in the clus- 
tering predictions with changes in cosmological parameters 
from their fiducial value. In this figure we varied the cos- 
mological parameters within their 2a limits determined by 
WMAP7 year data. The parameter in each case is fixed 
by fitting the observed luminosity function. All the curves in 
Fig. [S]are for z — 3 and threshold apparent magnitude 25. 
The parameter values which are different from the fiducial 
value are shown in each panel of Fig [5] The percentage dif- 
ference of clustering ranges from minimum of few % (for erg 
and Qb) and a maximum of ~ 12% (for n s ) for the assumed 
values of parameters. 

In particular we see that the correlation functions on 
large scales are not very sensitive to as and Qb- This is ex- 
pected for Qb as any small change in baryon density do not 
considerably alter the dark matter correlations, halo mass 
function and halo bias. On the other hand any change in erg 
will significantly modify all the above quantities. However 
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our models with various values of erg predict almost same 
large scales clustering. In order to understand why this ob- 
tains we consider two models A and B with different values 
of us given by ag = 0.80 and erf = 0.68. In Appendix A 
we have obtained an expression for the approximate ratio 
between galaxy correlation functions of models A and B on 
large scales (Eg. IA7[) . In particular we consider the cluster- 
ing of galaxies with threshold apparent magnitude 25. From 
our analysis in Section 0] we can obtain the average mass of 
a galaxy, M av , that can shine brighter than this magnitude 
threshold. For models A and B, the corresponding average 
masses are ~ 7.9 x lO 11 ^/© and M? v ~ 5.2 x 1O 11 M 
respectively. Using these in Eq. lA7l we get the ratio between 
correlation functions on large scales for models A and B to 
be 1.04. This is consistent with our model predictions shown 
in Fig. [5] Thus in our model various effects due to the change 
in erg cancel each other to get similar large scale clustering 
of high redshifts LBGs. 

The changes in all other cosmological parameters pro- 
duce notable differences in the large scale clustering espe- 
cially on very large scales. These modifications to clustering 
on large scales due to cosmology are quite interesting and 
promising. We have already seen that on large scales changes 
in astrophysical parameters do not strongly affect our model 
predictions of clustering as long as /* ft] is constrained by the 
observed LF. Moreover the percentage difference in correla- 
tion functions shown by various models in Fig. [S] appears to 
increase with the separation scale. Although the current ob- 
servations cannot distinguish between these models, future 
surveys with large number of sources detected in a large sur- 
vey volume may distinguish between large scale clustering 
predictions of high redshift LBGs with different cosmologi- 
cal parameters. Hence we conclude that our model for large 
scale correlation function of high redshift LBGs can be po- 
tentially used to constrain cosmological parameters. 



4 THE CORRELATION FUNCTIONS 
INCLUDING THE SATELLITES 

In the previous sections we assumed that a halo can host 
at most one detectable galaxy. This assumption is not ad- 
equate to explain clustering at small angular scales, espe- 
cially on scales smaller than virial radius. Each halo can 
in principle host multiple galaxies. A complete description 
of the distribution of galaxies inside a halo is called halo 
occupation dist r ibutio n (HOD). Following the approach of 
iKravtsov et alj (|2004T ) we separate t he central an d sate l- 
lite contribution s to HOD (see also [Zheng et al.l (|2005l ); 
ICoorav fc Ouchil (|2006l ); IConrov et al.1 (|2006h 1. That is the 
mean number of galaxies inside a dark matter halo can 
be written as N g (M) = f cen (M) + N S (M) where f cm {M) 
and N S (M) are respectively the mean number of central 
and satellite galaxies. We also assume the central galaxy to 
be situated at the center of the h alo and satellite galax- 
ies around it (|Kravtsov et alj|2004f ). In this approach the 
1-halo term has contributions from the correlation between 
central-satellite and satellite-satellite pairs. 

In the framework of this app roach, the the total cor re- 
lation function can be written as (|Coorav fc Shetbll2002T i 



where each term on RHS has contributions from central as 
well as satellite galaxies. 



4.1 The HOD: a physical model 

Often the HOD is modelled in a parametrized form, moti- 
vated by simulations or observations, and the value of the 
parameters are derived by fitting the galaxy correlation func- 
tion. Instead we adopt a more physical approach. We ask, 
given our model for galaxy formation, what is the expected 
mean occupation number of central and satellite galaxies of 
a given luminosity, in a given halo. These can then be used 
to compute the correlation function. 

We begin by formally expressing the mean occupation 
number of central galaxies with a luminosity threshold Lth 
inside a halo of mass M at any redshift z. This is given by, 



fcen(L t h,M,z) 



J dn(M,z c ) Q {L ^ ^ _ ^ ^ 
2 dz c 

°? dn(M, z c ) 



dz c 



z r 2 dn{M, Zc) 
dz c 



d.Z c 



n(M,z) 



(14) 



Ur) = e 3 h (r)+e g h (r) 



(13) 



Thus, in our model f C en(L t h, M, z) is the probability that a 
halo of mass M at z hosts a galaxy with brightness greater 
than Lth - Again the O (Lth(M, t(z) — t(z c )) function ensures 
that only those galaxies collapsing at z c and having L > 
L t h at z are counted in the integral over z c . As described 
earlier, zi(M, Lth, z) and zi(M, Lth, z) are the two redshifts 
at which a galaxy of mass M has to be formed, so that it 
shines with an observed luminosity Lth at z. 

In Fig. [6] we have plotted as thin lines, the average oc- 
cupation number of central galaxies (f C en(M)) as function 
of the mass of the parent halo, calculated using the above 
prescription. The results are shown at z — 3, 4 and 5 for 
three different magnitude thresholds which are labelled in 
each panel. These curves are obtained using the fiducial set 
of model parameters (given in Section 2) that reproduce the 
observed luminosity functiorQ. We can see that for a given 
redshift and limiting magnitude f C en{M) increases with the 
mass of the hosting halo to a limiting value of 1. Also note 
that the mean occupation number of central galaxies drops 
to zero below some mass scale M m m- This value of M m in is 
almost the same as that given in Table. [1] the change is due 
to the minor corrections to f,/rj we apply to take care of 
the satellite contribution to the LF. In Table [2] we give the 
values of Mmin for all the three redshifts and limiting mag- 
nitudes. For example, at z = 4 for apparent magnitudes 25, 
25.5 and 26, M min = 4.8 x 10 11 , 3.1 x 10 11 and 1.9 x 10 11 M© 
respectively. Many of the se mi analytical model s for cluster- 
ing o f low redshift galaxies (|Zehavi et al.l 1201 ll ; iMore et al.l 
l2009h use a step like function with smooth cut-off profile for 
fcen(M). One can see from Fig.[6]such a profile for f ce n{M) 



1 Inclusion of the contribution of satellite galaxies slightly mod- 
ifies the luminosity function and hence the best fit f*/rj (see be- 
low). In all the subsequent calculations we use this new best fit 
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naturally arises from our simple physical model of star for- 
mation. Infact the smoothness of the cut-off seen in the func- 
tion f CS n(M) at the low mass end for our models, is due to 
the scatter in the M-L relationship. We also give in Table [2l 
the average mass M av of a halo hosting the central LBG, for 
each redshift and luminosity threshold. This average mass 
is defined by, 



M a 



/ °°dM Mf cen (M,z)n(M,z) 
J™dMf cen (M,z)n(M,z) 



(15) 



4 and for apparent magnitudes 25, 25.5 and 26, 
ill c „ „ ln ii and 3 7 x 1Q ii Mq respectively. 

1.4 while at 



At z 

M av = 8.7 x 10 11 , 5.7 x 10 1 
At z = 5 these masses are smaller by a factor 
z = 3, M a v is larger by a similar factor. 

We also note that b g (k, Lth, z) in Eq. ([9]) can be rewrit- 
ten in terms of f C en(Lth, M, z) as, 

1 f x 

b g {k,L th ,z)= — r/ dM n{M,z)b{M,z) x 



n g (L t 

fcen(L th ,M,z)u(k,M,z) 



(16) 



which is the standard formula for the galaxy bias, excluding 
the satellite contribution. 

In order to get an estimate of the number of subhalos 
and thereby the number of satellite galaxies hosted by a 
halo of mass M, we use the conditional or progenitor mass 
function. The conditional mass function gives the comoving 
number density of subhalos of mass M s which formed at z c 
inside a region containing a mass M (or comoving volume 
V) that have a non - linear over density S a t z. It is given by 
i|Mo fc Whitdl 19961 ; ICoorav fc Shethll2002D 



n(M; 



,z c \M,S,z)dM s = f£v 10 f(„ 10 )^ 

M s V\o 



where 



Via 



(SdzJ-S^z)) 2 



! (M S 



\M) 



(17) 



(18) 



Here, Si(S, z) is the linear density contrast at redshift z 
corresponding to non-linear over density S and p m = 
QmPcric is the back ground density of baryons and cold 
dark matter. It should be noted that z c is greater than 
z and M s is smaller than M. The function on the RHS 
of Eq. (|17[) has the same form as the unconditional mass 
function but with 8 2 (z c )/o 2 (M) being replaced by (8 c (z c ) — 
6i(6, z)) 2 /{a 2 (M 3 )-a 2 (M)). We use Sheth-Tormen form for 
the RHS in Eq. (1171). in our calculations (see Section 3.3 of 
ICoorav fc Shethl 12002) for more details). The total num- 
ber halos of mass M s which formed at z c inside the volume 
containing mass M is 

N{M a ,z c \M,S,z)dMs = (M/p m )n{M s ,z c \M,5,z)dMs, 

(19) 

where we have multiplied n by the comoving volume 
(M/ Pm ). 

Our aim is find the number of satellite galaxies of a 
particular luminosity inside a dark matter halo of mass M. 
To calculate this we assume that the region containing mass 
M (Eq. (|17[) ) is a collapsed dark m atter halo at z. Thus i n 
Eq. dTTJ) we take 5i(5,z) = S c {z) (|Coorav fc Shethl |2002| ). 
In this limit the conditional mass function gives the number 
density of satellite halos in the mass range M s and M s + dM 3 
at z c inside a halo of mass M that collapsed at z. We choose 
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Figure 6. The halo occupation distribution {{N g (M))) as a 
function of the mass of the hosting halo, as predicted by our 
models at various rcdshifts for three limiting magnitudes. In each 
panel the thin curves corresponds to fcen(M) where as the thick 
curves give the total occupation N g (M) = f ce „(M)+N s (M). The 
redshifts and apparent magnitudes are labelled in each panel. The 
curves are obtained using a fiducial set of model parameters (given 
in Section 2) that reproduce the observed luminosity function. In 
addition to that, to obtain the satellite occupation (N S (M)), we 
have adopted Ato = 1.35 X td yn ( see text for details). 



the time derivative of N(M S , z c \M, S, z) as the formation 
rate of subhalos of mass M s inside a big halo of mass M. 
This can be justified by the the assumption that even though 
dark matter halos are formed and destroyed inside the over 
dense cell of mass M, the satellite galaxies formed in these 
mini halos have survived and can be observed. 

We further assume that no subhalo should be formed 
very close to the formation epoch of parent halo. Or more 
precisely if t c is the age of the universe when a parent halo 
collapsed then all the subhalos formed inside that parent 
halo within a time interval Ato prior to t c do not host a 
satellite galaxy. Rather they will be part of the parent halo 
itself. In our models we vary Ato as a free parameter. Also 
we expect that it is of the order of dynamical time scale 
tdyn of the parent halo. It should be noted that the time 
scale Ato corresponds to a redshift interval Azo which is a 
function of the redshift of collapse of the parent halcP]. 

2 We have also explored an alternate scheme whereby halos of 
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We will show that the typical mass of the halos that 
can host satellite galaxies satisfying the luminosity thresh- 
old condition is M > 10 12 M©. At high redshifts (2 > 3) 
most of these parent halos have collapsed close to the the 
redshift of observation. So for computational purpose we fix 
the redshift of collapse of big halos, that can host visible 
satellite galaxies, to be the redshift of observation (2) it- 
self. In the following discussions we also assume that star 
formation models for the central and satellite galaxies are 
identical (given by Eq. (0). This is similar to the assump- 
tion used by previou s authors (|Lee et al .11 20091 ; Ib errier et al.l 
120061 ; IConrov et al.|[200r3 ) where feedback due to halo merg- 
ing processes are ignored. 

We first compute the number of satellite galaxies of 
magnitude in the range Mab and Mab + cLMab inside a 
halo of mass M. This is given by, 
00 

dN(M s ,z c \M,z) 
<p(M,M AB ,z)dM AB = I v : 1 ' ' x 



dz 

z — Azo 

dMs dLi5oo 



dLl5Q0 dMAB 



dMAB dz c 



(20) 



It should be noted that ip(M, Mab , z)dM AB is not the lu- 
minosity function of satellite galaxies, but the total num- 
ber of satellite galaxies in a magnitude range Mab and 
Mab +dMAB inside a halo of mass M. The average number 
of satellites, N s , with luminosity greater than L t h in halo of 
mass M at redshift 2 can be computed as 



(N s \Lth, M, z)) = / <p(M,L,z)dL 



(21) 



where the luminosity L is related to Mab in the standard 
manner. 

We have plotted in Fig. [6] the total occupation num- 
ber of galaxies, N g (M) = f C en(M) + N B (M), as a function 
of the mass of the parent halo. These curves are shown in 
thick lines for at z = 3, 4 and 5 and for three luminosity 
thresholds. Thin lines correspond to occupation number of 
central galaxies. In order to obtain N S (M), apart from the 
fiducial model parameters that fit the observed luminosity 
function, we have adopted Ato = 1.35td yn - We find that for 
each limiting magnitude the mean number of satellites is a 
monotonically increasing function of hosting halo mass. As 
an example at z — 4 for a halo of mass 2 x 1O 12 M0 the mean 
satellite occupation numbers are respectively 0.20, 0.33 and 
0.44 at threshold apparent magnitudes of 25, 25.5 and 26. 
For a bigger halo of mass 5 x 10 12 Mq these numbers change 
to 0.43, 0.57 and 0.69 for the same limiting magnitudes at 
the same redshift. Thus the number of detectable satellites 
is typically less than unity. This implies that all halos do not 
host an additional detectable satellite galaxy; however some 
fraction of them do and it is these pairs of LBGs which con- 
tribute to the small scale clustering. Such a conclusion has 



satellite galaxies must have a mass smalle r than the parent halo by 
a factor f m (same as the parameter <p in lLee e~a l. (2009)). Since 
smaller mass halos typically collapse earlier than larger masses, 
the parameter Ato used here achieves a similar goal to f m ; al- 
though it does not exclude similar but smaller mass satellites in 
the main parent halo. 



also been arrived at bv lConrov et al.l (|2006l i using numerical 
simulations. 

It is also of interest to ask for a characteristic mass Mp 
of a parent halo which is hosting a detectable satelite galaxy. 
We define this as follows: 

_ J °°dM MN s (L t h, M, z)n(M, z) 
P f °° dMN s (L th ,M, z)n(M, z) ( ' 

This is also given in Table [2] For apparent magnitudes 25, 
25.5 and 26, and 2 = 4, M P = 1.6 x 1O 12 M , 1.0 x 10 12 M Q 
and and 6.3 x 10 11 M Q respectively. For z — 3 the corre- 
sponding Mp is ~ 1.5 times larger, while for 2 = 5, Mp is 
smaller by a factor ~ 1.4. This also means that the charac- 
teristic mass of parent halos hosting detectable satellites is 
~ 2 times larger than the mass of the halo hosting a central 
LBG. 

Note that the mean satellite number obtained in our 
physically motivated models can be fit asymptotically by a 
power law of the form N S (M) cx M a , with 01 ~ 0.7 - 0.9. 
Interestingly t his is very simi l ar to the parametrized form 
N S {M) used in lHamana et all (|2006h . who also adopt a sim- 
ilar slope. Moreover, the form of HOD that we derive from 
our physically motivated model ( and shown in Fig. O is also 
very similar to that obtained by IConrov et all (|2006h from 
their models of HOD using N-body simulations combined 
with a prescription of associating mass to light (see Figures 
8 and 10 of their paper). 



4.2 The one halo term 

We can now compute the 1-halo term using the standard as- 
sumption that radial distribution of satellite galaxies inside 
a parent halo follow t he dark matter density distribution 
(|Coorav fc Shethll2002l ). For the present calculations we use 
the NFW form for the dark matter de nsity distribu t ion. I n 
this case the 1-halo ter m is given by (iTinker et al.l d2005l) ; 
ICoorav fc Ouchl l|2006h : see also lSkibba fc Shethl (|2009l V) 



£ (L t h,R, z) 



I dM n(M,z) x \2f cen (Lth,M,z) 



(N s \L th ,M,z) pNFW ! Rm + 



(N s (N s -l)\L th ,M, 2) 



M 

\nfw(R\M) 
M 2 



}■ (23) 



Here 



n T g {L th ,z) = J dMn(M,z)(f cen (M,z) + (N s \L th ,M,z)) 

(24) 

is the total number density of galaxies which includes both 
the central and satellite galaxies. Further pnfw is the NFW 
profile of dark matter density inside a collapsed halo and 
^NF w( r \M) is the con volution of this density profile with 
itself l|Sheth et al.ll200lh . For the halo co ncentration param - 
eter we use the fitting functions given by (jPrada et al.l2012h . 
These fitting functions of concentration parameter at high 
redshifts, unlike earlier fits, flatten and then upturn with 
increasing mass. We also check the sensitivity of our results 
to other expressions for concentration parameter. Also fol- 
lowi ng the N-body simula t ions and semi-anal ytical models 
fee. iKravtsov et al.l i|2004l ), IZheng et all (120051 )). we assume 
that the number of satellites inside a parent halo forms a 
Poisson distribution. Thus we have {N S (N S - 1)) = TV 2 . 
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Table 2. Physical parameters constrained using observed lu- 
minosity function of high-z LBGs after incorporating the satel- 
lite galaxies in our semi-analytical models. These parameters are 
f*/Tj, an indicator of the light to mass ratio at any redshift, 
M m i n (Lt,hi z ), the minimum mass of a galaxy that can shine 
brighter than a given luminosity threshold L t ^ at redshift of ob- 
servation, M av (L t i l , z), the average mass of a dark matter halo 
that can host a central LBG satisfying the luminosity thresh- 
old Lthi ^pij^th, z ) an d the characteristic mass of a parent dark 
matter halo that hosts satellite LBGs which satisfy the luminosity 
threshold L t ^. 



z 




m 


M min /M e 


Mav/M® 


M P /M e 






24.5 


6.3 x 10 11 


1.2 x 10 12 


2.4 x 10 12 


3 


0.037 


25.0 


3.8 x 10 11 


7.9 x 10 11 


1.4 x 10 12 






25.5 


2.5 x 10 11 


5.2 x 10 11 


8.9 x 10 11 






25.0 


4.8 x 10 11 


8.7 x 10 11 


1.6 x 10 12 


■1 


0.035 


25.5 


3.1 x 10 11 


5.7 x 10 11 


1.0 x 10 12 






26.0 


1.9 x 10 11 


3.7 x 10 11 


6.3 x 10 11 






25.5 


4.3 x 10 11 


7.1 x 10 11 


1.3 x 10 12 


5 


0.032 


26.0 


2.7 x 10 11 


4.7 x 10 11 


8.3 x 10 11 






26.5 


1.7 x 10 11 


3.1 x 10 11 


5.3 x 10 11 



t — I 1 1 ■ 1 ■ r 




22 -20 -18 -16 

-m 

Figure 7. The best fit luminosity function of LBGs at z = 4 
(solid-black) including the contribution from satellite galaxies. 
Also shown is the individual luminosity functions due to central 
(dotted-red) as well as satellite (dashed-blue) galaxies. The ver- 
tical (solid-black) lines are drawn at absolute magnitudes -21.0 
and -21.6. 



4.3 The two halo term 

The correlation between satellite-satellite and satellite- 
central pairs located in different halos also modifies the 2- 
halo term, fn this case the expression for scale dependent 
galaxy bias in Eq. (|16p can be modified by adding in the 
contribution from satellite galaxies. Note that each halo of 
mass M has N 3 (Lth, M, z) satellite galaxies, brighter than 
L t h- Thus the number density of satellite galaxies associ- 
ated with a halo in any mass interval M and M + dM is 
n(M, z)N 3 (Lth, M, z)dM . These satellites are residing in a 
parent halo with a large scale bias b(M,z). Therefore, the 
average galaxy bias obtained by adding the contributions of 
both the central and satellite galaxies is, 

b g (k,z,L th ) = / dMn(M,z)b{M,z) 

n a {Lth, z) J 

[Un(L th , M, z) + (N s \L th , M, z))] u(k, M, z). (25) 

If halos do not host any satellite galaxies, ie 

(N s \L t h, M, z)) — 0, this expression reduces back to 

Eq. (|16l) . We can now compute the two halo term by 
substituting Eq. (f25]> into Eq. (g]) and using Eq. (fTTj) . 



In Fig. [7J we have plotted the luminosity functions of 
both central and satellite galaxies separately at z — 4. It is 
clear that the number density of satellite galaxies (dashed 
blue curve) at any given Mab is roughly an order of mag- 
nitude less compared to that of central galaxies (red dotted 
curve) at any luminosity bin. Thus the contribution to to- 
tal luminosity function (solid black curve) due to satellite 
galaxies is negligible at any given luminosity. However, these 
satellites can still contribute to the small scale clustering. As 
discussed earlier, the typical mass Mp of parent halos host- 
ing satellite galaxies is roughly 2 times the average halo mass 
Mav of a detectable LBG itself. As the luminosity roughly 
scales with halo mass, the central galaxy in the parent halo 
could then be about 2 times (or about 0.60 magnitudes) 
brighter. From Fig. [7] we can see that satellite galaxies with 
absolute magnitude say —21, in a parent halo of mass Mp, 
will be about half as abundant as the central galaxies, with 
absolute magnitude —21.6, in them. Such pairs which will 
then occur in every other parent halo of mass Mp will con- 
tribute significantly to the small angular scale clustering. 



4.4 The total luminosity function 

One may wonder if the additional satellite galaxies inside 
a halo can modify the luminosity function. In this section 
we compute the contribution to galaxy luminosity function 
from visible satellite galaxies. The total luminosity function 
of LBGs at any redshift can be written as, 

<f>{M A B,z) = ®cen(MAB,z) + $ sat (M AB , z). (26) 

In section [2] we computed the luminosity function due to 
central galaxies. The contribution to luminosity function due 
to satellite galaxies can be written as, 

$>sat{M A B,z)dM AB = f dM^<p(M,M AB ,z)n(M,z). 

(27) 
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5 COMPARISONS WITH OBSERVATIONS 

5.1 The total correlation functions 

The total two point correlation function of galaxies at any 
scale is the sum of the contributions from one halo and two 
halo terms (as given in Eq. {13])). The angular correlation 
functions can be estimated, using Eq. (|12p where one has to 
use this total galaxy two point correlation functions, as de- 
scribed in previous sub-sections. We again adopted the same 
fiducial model with n — 1.0 and three different M agn values, 
l.Ox 10 12 M©, 1.8x 10 12 A/ Q , and 2. Ox 10 12 A/ Q at redshifts 3, 
4 and 5 respectively. The cosmological parameters are kept 
to their fiducial value. We also use Ato = 1.35 td yn for all the 
models, but test the sensitivity of the results to this choice. 
This single free parameter will correspond to different red- 
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Figure 8. The angular correlation function of LBGs, taking into account of one halo and two halo terms, at redshifts 3, 4 and 5 for 
various limiting magnitudes. Each row corresponds to a particular redshift. which is labelled in the first panel of that row. In each 
row there are three panels showing the clustering predictions for galaxies with three threshold magnitudes, that are also labelled in the 
respective panels. Our fiducial model predictions of galaxy angular correlation functions with Ato = 1.35 td yn is shown in blue solid line. 
The dotted and dashed curves are for Ato = 1-0 td yn and 1.7 td yn respectively. The black dash-dotted lines show the co rrelation functions 
compu ted without satellite contribution (see Fig. [3]> - The data poin ts and error bars shown by black triangles are from lHildebrandt et alj 
(2009) and those shown by red squares are from lOuchi et alj 1 120051) . 



shift intervals Azo at z = 3, 4 and 5. Various parameters of 
our fiducial model are tabulated in Table. 

The total angular correlation functions computed us- 
ing our prescription for three redshifts and three threshold 
magnitudes are overplotted on the observed data in Fig. [8] 
In this figure the blue solid curves are our predictions of 
total galaxy angular correlation functions for the fiducial 
model with Aio = 1.35 tdyn- The dotted and dashed curves 
are for Ato = 1.0 td yn and 1.7 tdyn respectively. We see 
that adopting Ato = 1.35 tdyn provides a better fit to the 
data for all redshifts and magnitudes. A smaller (larger) Ato 
leads generically to an excess (deficit) of small angular scale 
clustering. The black dash-dotted lines show the contribu- 
tion to correlation function due to central galaxies alone 
(computed in Section [3] and given in Fig. |3J. The satellite 
galaxies contribute negligibly to clustering at large angu- 



lar scales, as suggested earlier in Section 3. The remaining 
contribution, which dominates especially at small angular 
scales (8 < 10"), is mainly due to the one halo term. It is 
also clear that the one halo term does not affect clustering 
on large angular scales (6 > 10"). Note that the one halo 
and two halo contributions are distinctly seen in the ob- 
served data for z — 3 and z — 4 as predicted by our models. 
But these distinct contributions are not so clearly seen in 
the z — 5 data. 

From the figure it is clear that our simple physical model 
gives a reasonable fit to the observed angular correlation 
functions at large (6 > 80") and small (8 < 10") angu- 
lar scales. This is true at all magnitude thresholds and at 
all redshifts. Moreover at z — 3 the the predicted angular 
correlation functions fits the observed data reasonably well 
at all angular scales (including the range 10" < 9 < 80"), 
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for magnitude thresholds 25.0 and 25.5. However, for the 
most luminous galaxies at z = 3 (with m < 24.5), and for 
redshifts 4 and 5, there is a discrepancy in the interme- 
diate scales between the theoretically predicted correlation 
function and observed data points. This could be due to 
(i) uncertainties related to distribution of satellite galaxies 
that is assumed to follow the NFW density profile of dark 
matter halos, (ii) the non-linear cor rections to dark matter 
power spectrum (|Mo fc Whitelll996f) and (in) the hig her or- 
der corrections to galaxy bias (|Coorav fc Shethll2002r ) which 
is not incorporated in our model. We discuss these possibil- 
ities further below. Nevertheless, it is quite remarkable that 
the predictions from our simple physical model fit the ob- 
served clustering over a wide range of scales and for the full 
range of z = 3 — 5. 




1 lO lOO lOOO 



5.2 Redshift evolution of the clustering 

In the previous section we focused on the clustering at red- 
shifts 3, 4 and 5. In order to see the clustering evolution of 
LBGs we have plotted in Fig. [9] our model predictions of 
angular correlation functions at redshifts 3, 4 and 5 for an 
absolute magnitude ~ —20.1 (corresponding apparent mag- 
nitudes at these three redshifts are 25.5, 26 and 26.5 respec- 
tively) . From the data points it is clear that at on large scales 
(6 ^ 80") the clustering is maximum at z=5, decreases con- 
siderably by z = 4 and is almost the same at z = 3. Interest- 
ingly one can see from the figure that our model predictions 
of correlation functions also show such a trend. Note that 
the magnitude of the statistical bias is expected to decrease 
with decreasing redshift, for a fixed halo mass. On the other 
hand the dark matter power spectrum will increase with de- 
creasing redshift. Thus the clustering amplitude of halos is 
expected to first decrease with z and then subsequently in- 
crease. How exactly galaxy clustering evolves with z depends 
also on the mapping between the mass of a dark matter halo 
and the luminosity of the galaxy it hosts. It is interesting 
that our simple physical model gives a mapping between M 
and L which qualitatively predicts the observed redshift evo- 
lution of clustering on large angular scales. On small scales, 
where the one halo term is dominant, our fiducial model 
is also able to qualitatively reproduce the observed redshift 
evolution of the angular correlation functions. 



5.3 Sensitivity to the halo density profile 

The discrepancy noted above, between the observed and our 
model predicted w(0) at intermediate angular scales, could 
be due to one halo correlations of galaxies extending beyond 
the virial diameter of NFW density profile. Physically this 
means that subhalos can reside outside the virial diameter 
of the parent halo. If this is true, data suggest that, this 
effect is stronger for z ^ 5. 

To explore this possibility we assumed that halos follow 
a new density profile which is of the same form as the NFW 
profile, but with a larger virial radius f vir = s v i r r v i r . We 
keep the the concentration parameter and the mass of the 
parent halo the same as before. We took the fiducial values 
of s V ir = 1.5 and s V i r = 2.0 at redshifts 4 and 5 respectively 
(see Table [3}- Adopting s V i r > 1 can increase the clustering 
at intermediate scales but will reduce the amplitude at small 



Figure 9. The redshift evolution of clustering of LBGs for an 
absolute limiting magnitude -20.1. The data points and error bars 
are from Iffildcbrand t et al. We used the fiducial cosmol- 

ogy parameters in this figure. All the astrophysical parameters 
the same as in Fig. [8] 



Table 3. A summary of the fiducial parameters used in various 
models at at z = 4, 5 and 6. AGN feedback determined by the 
parameter M agn evolves with redshift even in our fiducial models. 
The quantity n that determines the duration of star formation in 
our models is fixed to be 1. Also we kept Ato to be a constant 
in the fiducial models. However for models where we extend one 
halo term beyond the virial diameter of a halo, we changed Ato 
slightly from fiducial value. 



2 


M agn /M @ 


Ato/tdyn 


Az 


Svir 


3 


1.0 x 10 12 


1.35 


0.80 


1.0 


i 


1.8 x 10 12 


1.35 


1.00 


1.0 




1.8 x 10 12 


1.25 


0.95 


1.5 


5 


2.0 x 10 12 


1.35 


1.25 


1.0 




2.0 x 10 12 


1.25 


1.20 


2.0 



scales. This can be compensated by adopting a smaller Ato. 
In Fig. [TO] we have shown the effect of extending the one 
halo term beyond the virial diameter of the parent halo. 
The value of Ato is also reduced to 1.25td yn to obtain a 
better fit at the smallest angular scales. The dotted black 
curve is for s„; r = 1.0 (fiducial model). The curves obtained 
after extending the one halo term beyond the virial diameter 
of dark matter halos are shown in solid red lines. One can 
clearly see that these new curves provide a better fit to the 
observed clustering in small angular scales (6 < 20"). It 
is also clear from the figure that a change in the form of 
distribution of satellite galaxies inside the halos alone will 
not explain the discrepancy seen at 20" 9 ^ 80". 

We have also considered the prescription for the con- 
centration parameter given bv lBullock et al.l (l200ll). which is 
used by earlier works high redshift clustering ( Hamana et alj 
I2006T) . The dashed blue curves in Fig. [10] show the cor- 
responding results at z = 4 and 5. One can clearly 
see that adopting the concentration parameter given by 
iBullock et al.l l|200ll 'l underpredicts the clustering at 6 ^ 3" 



© 0000 RAS, MNRAS 000, 000-000 



14 C. Jose et al. 



o. i — 



O.Ol 



10 



O. 1 



O-Ol — 




•3- ( 1 I t Ht#f H iHWHf 



1 O 1 OO 1 O 3 1 



io lOO 
©(sec) 



IO 3 1 



IO 



lOO 



IO 3 



Figure 10. Effect of clustering predictions on concentration parameter and s„; r Ou r fiducial model is sh own in black dotted curve. The 
dashed blue curve is for a model where one uses concentration parameter given by iBullock et alj feOOll) . The solid red line shows our 
model predictions, that show the effect of extending the one halo term beyond the virial radius of the parent halo. In these cases, for 
simplicity, we multiplied the virial radius by factor s v i r . For z = 5 we have used s v i r = 2 and for 2 = 4 we have used s v i r = 1.5. Note 
that our fiducial model corresponds to s„; r ~ 1. The data points and error bars are from lHildebrandt et al] j200gT) . We used the fiducial 
cosmology parameters in this figure. All the astrophysical parameters the same as in Fig. [8] 



5.4 Dependence on astrophysical and 
cosmological parameters 

In Section r3.ll we found that the large angular scale correla- 
tion functions (with 80" < 8 < 1300") were fairly insensitive 
to changes in the assumed astrophysical and cosmological 
parameters, provided these models still fit the observed LF. 
We now ask if this is true also for small angular scales be- 
low 8 < 10" where our fiducial models can reasonably fit 
the observed data. On small angular scales, the correlation 
functions can in principle be a function of astrophysical pa- 
rameters as these change the total number of central and 
satellite galaxies satisfying the given luminosity criteria. 

In Fig. 1111 we show the effect of varying k on small scale 
clustering. The blue solid curve is for our fiducial model with 
k = 1.0. The blue dotted and dashed curves are for k values 
0.7 and 1.3 respectively. All other parameters except for 
have been kept the same, and f,/rj is varied such that the 
predicted LF still fits the observed LF data reasonably well. 
One can clearly see that clustering predictions especially for 
k = 0.7 do not fit the correlation functions on small angular 
scales. Therefore we conclude that small scale clustering of 
LBGs at high redshifts can give useful information about k 
and thus the duration star formation in a halo. As a compar- 
ison, we note that, varying k in our model is like varying th e 
duty cycle in the semi-analytic models of iLee et alj l|2009h . 
We find that k ~ 1 or the star formation duration of order 
tdyn is favored by the small scale clustering data. 



We also varied M agn around the fiducial value at each 
redshift, ensuring that the predicted LFs still reasonably fits 
the and observed data. We find that even the small scale 
angular correlation functions are fairly insensitive to these 
changes. For example when we varied M agn at redshift 4 
from 1.5 x 10 12 Af Q to 2.1 x 10 12 Af©, the clustering predic- 
tions change at all scales by about —7% to 5% from the 
fiducial model. We also found that the situation is similar 
for all magnitudes and other redshifts. 

We have also considered the effect of varying a whole 
suite of cosmological parameters within their 2a limits de- 
termined by WMAP7 year data. Again we ensure that the 
predicted LFs best fits the observed data, by varying f*/f). 
Fig. \T%\ shows the percentage differences in the clustering 
predictions with changes in cosmological parameters from 
their fiducial value. All the curves are at z — 4 for a thresh- 
old apparent magnitude 25. 

From Fig. [12] one can see that the small angular scale 
clustering is very insensitive to changes in all the cosmolog- 
ical parameters except as. Infact the percentage difference 
in small scale clustering differs at most by ±4%. However 
changing as significantly alters small scale angular correla- 
tion functions. From the figure it is clear that the percentage 
differences in clustering predictions are in the range ±30% 
when as is varied between its 2a limits (0.68 to 0.92) given 
by WMAP 7 year data. Moreover the change is scale depen- 
dent, and thus varying as will change the shape of the small 
angular scale correlation function. 
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Figure 11. The change in angular correlation function with k around its fiducial value at z = 4. The blue solid curve is the fiducial 
model prediction with k = 1.0. The dotted and dashed curves are for k = 0.7 and 1.3 respectively. 





;.. " " " " r x " 

f ~- 


o, - 0.08 




i ■ 1 








n.'-0.*4 




i 


n. -O.S01 " 


1 o,h«- 6.082...". 






Q,h«- 0.023 




1 - 


h - 0.86 . 






h - 0.77 . 




1 ■- 


Q.h«-0.1063 






ft.h«-0.118fi 



10 



Figure 12. The percentage change is total angular correlation 
(small scale) function when one varies the cosmological parame- 
ters around the fiducial model at z = 4. Here the plots are for 
threshold apparent magnitude 25. The fiducial model is shown in 
blue solid curve. 



6 DISCUSSION AND CONCLUSIONS 

We have presented here a physically motivated semi- 
analytical model of galaxies to understand the clustering 
of high redshift Lyman break galaxies, where the model pa- 
rameters are constrained by the observed high z luminosity 
function. For this purpose we use and expand upon the stan- 
dard halo model. Galaxies are assumed to be formed inside 
dark matter halos. Their luminosity is determined by a phys- 
ical model of star formation, which is a function of the mass 
and age of the hosting halo. 

We began by assuming that each halo can host at most 



one visible galaxy. On fitting the observed LF, we determine 
the relationship between the luminosity of a galaxy and the 
mass of its host halo. This allows us to calculate the large 
scale bias for LBGs satisfying any luminosity threshold. This 
bias is then folded in with the dark matter power spectrum 
to predict the two point spatial correlation function of these 
LBGs and hence the angular correlation function w(0). For 
our fiducial model, we find that th e predicted w(9) compares 
very well with the observed data of lHildebrandt et alj (|2009h 
at 6 > 80", for the whole range of redshifts z = 3 — 5 and 
limiting luminosities (magnitudes) (see Fig. [3} . 

The predicted large scale galaxy b ias b g agrees well 
with t hat observationally determined by iHildebrandt et al.l 
|2009l ) (see Table [2). At a given z the bias increases with in- 
creasing luminosity, while for a given Lth it increases with z. 
However, we find a smaller sprea d in b q as a function of L t 



at any given z compared to that of lHildebrandt et alj (|2009t ) . 
Remarkably, we find that the predicted large scale cluster- 
ing of LBGs is fairly insensitive to the assumed astrophysical 
or cosmological parameters, provided we simultaneously fit 
the observed LF. This may point to an important internal 
consistency of our physical model; that if we fix the mass 
to light ratio correctly by using the LFs of LBGs, then the 
standard LCDM model correctly predicts the amplitude of 
their large scale clustering. 

We then extended our approach by incorporating the 
Halo occupation distribution, which provides the distribu- 
tion of galaxies inside dark matter halos. This is separated 
into the central, fcen, and satellite, N s , contributions. Of- 
ten the HOD is modelled in a parametrized form. Instead we 
have adopted a more physical approach. We use our prescrip- 
tion for computing the LF to estimate f cen . The conditional 
mass function and our star formation model are used to cal- 
culate the mean number and luminosity of satellite galax- 
ies in a parent halo. The calculated forms of / ce „(M, Lth, z) 
and N S (M, Lth, z) comp are reasonably with th at assumed in 
parametriz ed models of Hamana et al.l (|2006t ) and the sim- 
(|2006l ). At z = 4 and for apparent 



ulations of lConrov et al 



magnitude thresholds in the range 25-26, the average mass 
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of halos contributing to the observed clustering, M av , ranges 
from 8.7 - 3.7 x 10 n M o . At z = 5 these masses are smaller 
by a factor ~ 1.4, while at z — 3 these masses are larger by a 
similar factor. At any given redshift and magnitude thresh- 
old, the typical mass Mp of parent halos hosting detectable 
satellite galaxies, are about 2 times larger than M av . 

Having obtained the HOD, we can calculate both the 
one halo and two halo contributions to the total correlation 
function of LBGs. Our simple physical model gives a rea- 
sonable fit to the observed clustering of LBGs at z = 3 at 
all angular scales for magnitude thresholds 25.0 and 25.5. 
At other redshifts (z = 4 and z = 5), and for the most lumi- 
nous galaxies at z = 3 (with m < 24.5), the predicted w(9) 
again fits the observed data well at both large (9 > 80") 
and small (9 < 10") angular scales. The clustering at small 
angular scales is likely to be dominated by pairs of LBGs 
rather than rich clusters, as the number of detectable galax- 
ies hosted by most collapsed halos is typically less than 2. 

The small angular scale clustering is also not very sen- 
sitive to changes of several cosmological and astrophysical 
parameters from their fiducial values, as long as we simulta- 
neously fit the observed LF. However, the amplitude of w(9) 
on small angular scales is very sensitive to the value of k. 
There is also an increased sensitivity of both the amplitude 
and shape of w(9) at small angular scales to changes in ag. 
However, as ag is much better determined from other ob- 
servations, one may be able to constrain the duration of the 
star formation activity in LBGs using w(9) at small angular 
scales. The present data are consistent with n ~ 1 or a star 
formation duration of the order of the dynamical time scale 
tdyn of the dark matter halo. 

The theoretically predicted w{9) at intermediate angu- 
lar scales is smaller than that observed, and the reason for 
this remains to be fully understood. We explored in detail 
whether this excess can be due to a more extended satel- 
lite galaxy distribution. This only partly accounts for the 
discrepancy. We also find that the non-linear corrections to 
dark matter power spectrum does not raise the predicted 
clustering at the intermediate scales to the observed values. 
We suspect that the higher order corrections to galaxy bias 
could be playing a role, a possibility which we hope to ex- 
plore in the future. Nevertheless, it is noteworthy that the 
predictions from our simple physical model, employing only 
a few free parameters, can fit the observed clustering data 
over a wide range of scales, redshifts and limiting luminosi- 
ties. 
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bias is scale independent. Thus using Eq. (|A2|) we get 
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b g (z, L th ) 



M av 



(A3) 



The large scale (small k) dark matter power spectrum 
Pn n (k) always scale as square of erg. Hence using Eq. (|10|l 
and Eq. (|A3[l we get the follwing scaling law of galaxy power 
spectrum. 

2 ' 3 + "ef/> 

M 3 

P 2h {k,L th ,z) = b 2 g (Lth,z)Pn n (k,z) oc — ^—^ F(k,z) 

(A4) 

where Pi in (k,z) = a 2 F(k,z). The two point correlation 
functions on large scales also scale in this way. Thus under 
the above approximations 



M 3 

U(Lth,R,z)<x —22-= F(R,z) (A5) 



where F(R, z) is the transformation of F(k, z) as given in 
Eq. ()11[) Suppose one assumes two models A and B with 
different values of as, say a£ and a B . The ratio between two 
point correlation functions of both models on large scales 
will be 



f A 



i. 



Mg,, 



2 < 3 +"e//> 



(A6) 



For galaxies we have n e // ~ —2. Hence the ratio between 
galaxy two point correlation functions of models A and B is 
given by 



<i2h 
S2/! 



Mi 

Ml 



(A7) 



APPENDIX A: CLUSTERING DEPENDENCE 

ON a 8 

For any mass scale M, the variance of smoothed density 
contrast a 2 (M) oc a$k%[P(KM) ~ ogfc^~" e// . Here n e ff is 
the effective spectral index, which is ~ —2 on galatic scales 
and —1 on cluster scales. We also have k~j^ ~ M 1 ^ 3 . Thus 
we get 

a 2 {M) oc at M <3+ 3 " c "\ ( Al ) 

For high redshifts v(M,z) = (S c (z)/a(M)) 2 » 1. Hence us- 
ing Eq.|6]we get the halo bias, b(M, z) oc u(M, z) oc l/a 2 (M) 
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